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SUMMARY 

An analytical model for separated airfoil flows is presented which is 
based on experimentally observed physical phenomena. These include a free 
stagnation point aft of the airfoil and a standing vortex in the separated 
region. A computer program is described which iteratively matches the outer 
potential flow, the airfoil turbulent boundary layer, the separated jet en- 
trainment, mass conservation in the separated bubble, and the rear stagnation 
pressure. Separation location and pressure are not specified a priori. Re- 
sults are presented for surface pressure coefficient and compared with experi- 
ment for three angles of attack for a 6A(W)-1, 17% thick airfoil. 


INTRODUCTION 


Separated flow on wings has drawn researchers' interest ever since it was 
found to be responsible for aircraft stall. Interest has recently intensified 
due to the ability of the new GA(W) series i)f airfoils to operate stably with 
flow separated up to half of the upper surface length. However, research was 
long restricted to experimental observation due to the complexity of the phe- 
nomenon. For very low Reynolds numbers,, solution of the Navier-Stokes equa- 
tions should provide an analysis of the flow, but even with modern computers, 
this is a costly task. For realistic airfoils, turbulent flows occur and re- 
quire some sort of mathematical model. Several models have been proposed and 
solved by computers with various simplifying assumptions or empirical rela- 
tions. Most of the previously proposed models consider the separated region 
to extend to infinity downstream or it is modeled as a bulbous region. The 
actual physical processes in the separated region have not been included in 
the models. Although some of these methods have produced reasonably good com- 
putational results, their ranges of applications are generally limited and 
empirical input is usually required. 

A mathematical model was sought which included the principal physical fea- 
tures. The model was first formulated by logic, found to conform to the veloc- 
ity and pressure data, and then tested by special wind tunnel and electric ana- 

*This research was supported by NASA Grant NSG 1192 and is available in 
more complete form in NASA CR-145249 (ref. 1). « 
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log experiments. Finally, a computational program *Thif pro- 

tive of minimizing the empirical data were 

gram was used to predict pressures on an which detailed data were 

available for comparison and evaluation of the model. 

Symbols are defined in an appendix. 

ORIGIN OF THE FLOW MODEL 

A sketch of the flow model is shown in figure 1. The flow has separated 
from the upper surface and leaves the lower surface at the trailing edge. The 
let mixing^sheets starting from the two airfoil separation points coalesce to 
form thrieparulon bgbbll behind the airfoil. These two jets entrain air from 
trdeSS TeSion (near-wake). The entrained air 

flow of air which must be supplied from somewhere. If S and S (fig. 1) are 
tirtSo separation streamlines, it can be seen that the amount “i«J / 
flow entrained are growing in the <*<>""**Team dirMtion, consegire^^ 
cDace available for backflow decreases and the demand for it increases. :>ince 
Sfs%rnc}'ciStrnue very far. ^ tenni nation is r^ui red of^J wa e re- 

circulation region, and a stagnation point is formed. At this point tne two 
streams rejoin defining the end of the separation bubble. 

The entrained masses of the two jets aro not same, since their^l^^^^ 
and velocities are different. Generally, the upper one will entrain a lan^en 
mass. Therefore, the stagnating streamlines are not necessarily 
• "s and Fiaure 1 shows two oth*^r streamlines R and R stag 

iatina at ih^rw^ pStnt ^nd p^SJtding passages for the flow to enter and leave 
th^sIpaJation bubble. The mass flow >^ate through the 

s must be the same as that between R' and S'. Further, R and R do not ter 
minke at the stagnation point but must continue upstream and downstream. Mass 
rnservation in the regioS requires the formation of two vortices and an S- 
shaped corridor of flow from bottom to top. 

At this point experimental verification was sought for this seemingly log- 
ical model nS reference could be found to a stable free vortex above the 
trailing edge, nor was the existence of a free stagnation , 

At W Chita Stke University, wind tunnel tests were in process for the GA(W)-1 
airfoil Jrder NASA contract. Detailed measur^ents of P*^®J5“r® 
were made downstream of the trailing edge at high angles of attack 2) 

Fiqures 2 and 3 are typical results of these tests and can be seen ^ 

locally high pfessure region just aft of the trailing ^ge and (2) the 
recirculating flow region terminating at the rear high pressure point. 

Based on experimental data, constant pressure is assumed in the whole 
<;pnarated reqion except in the neighborhood of the stagnation point, where a 
locally high pressure region will form along the RR' line, turning bac e 
low velocity flows. 

Viscous effects can be ignored in the neighborhood of the rear stagnation 
point and all velocity changes considered as being due to the pressure gra 
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bo^®5l and pressure domi- 

boundary layers and more recentlv anni P»*andtl for 

attaching flows in the Chapman-K^rst mixin^modll^^^^ separated and re- 

poiit^ tCthr^ pressure region 

point. Th.s Plane is also the location 

determine its existenc^* A^^S^^crjhord’GAfW'!!'?^^ devised to 

speed tunnel with a thin aluminum "splitter" Placed in a low 

flow at midspan. The olate pytonHoH w« * mounted parallel to the 

length. The plate was pain tefo? spot °"® ®»’ord 

sta9nat'oS'Soint'’|Iju c*losrtJ°Jlle‘?rI?™„g edSe?'*" ^ 

vortex abSJe^thS’relrpoJtiS* of*the“ilrfSlK ° 

Join the main streak This S-sL jJ'* “I""'''' ‘o 

clearly seen on thrstreakld p?a' flow was 

separ. ted airfoil*fi^was''set up'"*'^Tto*eIect "P"*®'; electric analog of the 
adjusted to meet the Kutta condition P»‘P"‘’'o 1 of the airfiii was 

tional voltage point was aZe and «"■< «" n«<li- 

ated the S-shaped 'lower-tS-Spper"nw I'll ®dge. This cre- 

nation. The streamlines around the air?nn^ separation, and a rear stag- 
to the strength of the rea? potential So?nJ sensitive 

exact location. This fact was used inM,a^??„’“jS| 

COMPUTATION OF THE MODEL 

the airfo’h surface! ^nd^^ separated f[ow ^ boundary laye^- on 

etched, together with 

Assumptions Made 

1. Steady, incompressible, plane flow of air. 

2. The boundary layer is fully turbulent. 
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3. The rear stagnation point is located downstream of the trailing edge one- 
third the distance between the upper separation point and the trailing 
edge. This was empirically derived and the solution was found to be quite 
insensitive to this valuei as the electrical analog study suggested. 

4. There is a constant-velocity reverse flow region, a "core flow" for the 
recirculating air, having a velocity 20% that of the adjacent free stream. 
Experimental evidence for the constant-velocity profile was abundant but 
no logical model can be suggested to give a "core flow" from the reversal 
of a shear flow. Therefore, a purely empirical choice was accepted as 
necessary =*nd the 20% value derived from the examination of several 6A(W) 
wing flow measurements. Results were somewhat sensitive to this choice, 
and this is regarded as the one genuinely empirical value in the flow 
model computation. 

5. The streamlines which stagnate behind the body do so without loss of total 
pressure. This is the commonly used Chapman-Korst jet-mixing flow model. 
The pressure variation from the trailing edge to the stagnation point is 
approximated by a parabolic variation to conform to experimental trends. 
Results were found to be insensitive to this choice. 


Potential Flow 

A potential flow solution method was required which would give velocities 
and pressures for specified geometries on the forward and lower'airfoil sur- 
faces, and also produce streamline geometry for a specified pressure distribu- 
tion in the separated regions. The Mixed Boundary Condition potential flow 
program of McDonnell Aircraft Company, St. Louis, met this requirement and 
was made available. This program is a modification of the wing-body code de- 
veloped by Woodward (ref. 3). The airfoil is divided into a number of panels 
on the chordline and the thickness, camber and angle of attack are represented 
by planar sources and vortex singularities. Panel lengths chosen for computa- 
tion were 1% of chord near the leading and trailing edges, and these increased 
5% in the center of the airfoil. Near wake panels were also 1% long. The po- 
tential program treated the separated region as an extension of the airfoil, 
with the body terminating ct the rear stagnation point. 


Boundary Layer Analysis 

Head's entrainment method (ref. 4) of calcilating the turbulent boundary 
layer was chosen as being sufficiently accurate without undue complexity. The 
flow was assumed to be turbulent from the It-Jing edge or, on the lower side, 
from the front sfagnati^'n point. After the potential flow solution was accom- 
plished, the displacement thickness was computed and added to the airfoil sur- 
face. The augmented airfoil was again solved by the potential program and the 
process repeated to convergence. Since the potential program provided the 
,'eparation st.'eamline shape, the displacement thickness at the separation point 
was added to the ordinates of the separation streamline to obtain the effective 
displacement surf* re. 
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Tte crUeria for sepiratlon of the hoandarv layer ms the <aiM «r m.. 
shape factor, H. uhlch Is the ratio of the dlsplaceheot^ictae^s ^ 

™a»ntM iJi1ckn«s. Host Integral xethods assMe H 'aloS’SfTe to ?4 as 
the se^rat^on Indicator. The exact value appears to depend on surface cirva. 
Tiin* strer.^th, and possibly the upstream boundary layer h1stw*y 

15 L J!” “ • ■“» of ■omerlio' tte Sutton 

point to no,e. An Initial -.aloe of 2.? ms used »lth separJelSTarStJ^telld. 


Jet Nixing Analysis 

• ^ iniXinQ C0H^)Utdt1OI1 WftS fOf* the IIDD#ir ^^rnunr nnu^i 1 1 

fecJrof^Jhe''li;n?al°b'^" "«thod''of"llUrlf?*6ri« wed VS!idl^*thS^If" 


Determination of the Stagnating Streamlines 

ouirJi^ts^^.rrSfr'!!:!?/ *?“ ?;.<»« «9- «) «" *« fro« Ue t«o re- 

Ke sepa Stioi TolSi M S!5t‘ STueTISerr"*’ 

leaving between streamlines S and R. The model assumes a constant statiJnroc 

KM 

located . ^ recalculated and the new R and R* Streamllies 


Recirculating Mass Balance 


.oca/55Z?,‘,%?5; itrx ?rs:‘?e;™n’cVi tz :5Joc’fS“’?hVMr 
Sru‘5[<i; 55"^ ~ ??* 
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COMPUTED RESULTS 


Figure 7 lllustretes the program logic. In this figure, Pp refers to the 
pressure at the rear sUgnatlon point. This program <^s us^ to calculate the 
pressure distributions on the 6A(ii)-l 17* thick airfoil at three «" 9 les of at- 
tack (18.4*. 16.4* , and 14.4*) for which experimental data were available 
(rff. 7). The results are shown in figures 8. 9 and 10. It can ^ seen that 
agreement with experiment is good. The separation pressure is predicted quite 
aLurately and the separation location is slightly aft of ihe 
values. The rear stagnation pressures also agree well with experiments, com- 
parative values are given in Table 1. 

The method does not prescribe any of the separaiicn variables. Examina- 
tion of the computational histories shows that the sewration ^int and pres- 
sure vary freely as the iterations proceed. Accuracy is dependent upon the 
number of panels used. The computed pressures show deviations from ex^ri- 
mental values at about 15* chord where panel size was Increased suddenl^y. By 
choosing a larger number of panels the accuracy would ^ improved at the cost 
of increased computing time. Using 29 panels on the airfoil, computation times 
on an IBM 360/44 or 370/145 were 14 to 40 minutes, depending on angle of at- 
tack and iterative accuracy desired. 

Some improvements are needed in the computational elements. The 
program used represented the airfoil by a singularity distribution on the chord 
line. This is adequate for thin airfoils, but for airfoils with large thick- 
ness or angle of Stack, there can be significant errors. A potential pi^ram 
using surface singularities (but still of the mixed boundary layer condition 
type) would improve accuracy. Similarly, Head's boundary layer method is us- 
able and simple, but is not the most accurate available. Computations should 
be made for other airfoils to evaluate the reverse flow velocity assimption. 


However, even at this point of development, this model has been shown to 
include all of the significant physical features, to be ca,.^Jle of computation 
with reasonable computer times, and to give good surface pressure results. 


f 
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APPENDIX 


SYHBOLS 



H 

M 

M 

Pr 

R.N. 

R 

S 

^sep 


Chord Length 
Pressure Coefficient 

Pressure Coefficient at Rear Stagnation Point 

Shape Factor of Boundary Layer 
Mach Number 
Mass Flow Rate 

Pressure at Rear Stagnation Point 

Reynolds Number Based on Airfoil Chord 

Streamline which stagnates at rear 

Streamline wh^ch separates from surface 

Length from separation point to airfoil trail inq 
edge; fraction of chord length. 


a Angle of Attack 

5 Boundary Layer Thickness 

6 * Boundary Layer Displacement Thickness 

o Jet Spreading Parameter 


Subscripts : 

L Lower side of airfoil 

SEP Separation point 

U Upper side of airfoil 

“> Free stream value 


Superscript ; 

Prime indicates lower-side flow 
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Figure 1.- Schematic diagram of separated flow model. 


ALPHA =18 -4 DEG- 
RN = 0-82EE 07 
tMCH ND- 0-130 



Figure 2.- Experimental velocity plot. GA(W)-1 airfoil. 



Figure 3.- Static-pressure field contours. GA(W)-1 airfoil 
a - 18.4"; R.N. « 2.2 x iq 6 (from ref. 2). 
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Figure 7.- Computer 

program logic. 





















